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o 

, We study the directional dichroic absorption spectra in the x-ray region in a polar ferri- 

magnet GaFeOa. The directional dichroism on the absorption spectra at the Fe pre-ii'-edge 
, arises from the EI-E2 interference process through the hybridization between the Ap and 

3d states in the noncentrosymmetric environment of Fe atoms. We perform a microscopic 
calculation of the spectra on a model of FeOg with reasonable parameter values for Coulomb 
interaction and hybridizations. We obtain the difference in the absorption coefficients when 
the magnetic field is applied parallel and antiparallel to the c axis. The spectra thus ob- 
^ , tained have similar shapes to the experimental curves as a function of photon energy in the 

^ , Fe pre-if-edge region, although they have opposite signs. 

KEYWORDS: directional dichroism, magnetoelectric effect, GaFeOs, noncentrosymmetricity, 

' O ' E1-E2 interference, x-ray K-edge absorption 
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1. Introduction 

Gallium ferrate (GaFeOs) exhibits simultaneously spontaneous electric polarization and 
' magnetization at low temperatures. This compound was first synthesized by Remeika,^^ and 

C5 I the large magnetoelectric effect was observed by Rado.^^ Recently, untwinned large single 

, crystals have been prepared, and the optical and x-ray absorption measurements have been 

■ carried out with changing the direction of magnetization.^'^) The purpose of this paper is 

to analyze the magnetoelectric effects on the x-ray absorption spectra and to elucidate the 
microscopic origin by carrying out a microscopic calculation of the spectra. 

The crystal of GaFeOs has an orthorhombic unit cell with the space group Pc2in.^^ 
Each Fe atom is octahedrally surrounded by O atoms. With neglecting slight distortions of 
octahedrons, Fe atoms are regarded as slightly displaced from the center of the octahedron; the 
shift is 0.26A at Fel sites and — O.llA at Fe2 sites along the b axis.'^^ Thereby the spontaneous 
electric polarization is generated along the b axis. Note that two kinds of FeOg clusters exist 
for both Fel and Fe2 sites, one of which is given by rotating the other by an angle vr around the 
b axis. As regards magnetic properties, the compound behaves as a ferrimagnet,^'^) with the 
local magnetic moments at Fel and Fe2 sites aligning antiferromagnetically along the c axis. 
One reason for the ferrimagnetism may be that the Fe occupation at Fel and Fe2 sites are 
slightly different from each other. ^) In the present analysis, we neglect such a small deviation 
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photon wave vector 



M loc : local magnetization 




Fig. 1. Geometry of absorption experiment."''^ X-rays propagate along the a axis with polarization 
along the b axis or the c axis. The electrie clipolc moment is along the b axis. When the magnetic 
field is applied to the positive direction of the c axis, the sublattice magnetization is directed to 
the negative and positive directions of the c axis at Fel sites and at Fe2 sites, respectively. When 
the magnetic field is reversed, the sublattice magnetization is reversed. 

from a perfect antiferromagnet. 

In the K-edge absorption experiment,^) the x-rays propagated along the positive direction 
of the a axis, and the magnetic field was applied along the ibc axis, as illustrated in Fig. 1. The 
difference of the absorption coefficient was measured between the two directions of magnetic 
field. As will be shown later [eq. (4.11)], these difference spectra have characteristic dependence 
on the polarization and magnetization, and would be termed as magnetoelectric spectra.^) 
Since the compound is a ferrimanget, reversing the direction of applied magnetic field results 
in reversing the direction of the local magnetic moment of Fe atoms. 

In the analysis of absorption spectra, we consider only the processes on Fe atoms, since the 
Is-core state is well localized on Fe atoms. In addition to the El-El and E2-E2 processes, we 
formulate the E1-E2 interference process which gives rise to the magnetoelectric spectra. The 
contribution of this process arises from the mixing of the Sc/^ls-configuration to the ApSdPls- 
configuration, where Is indicates the presence of a Is-core hole. Such mixings exist only under 
the noncentrosymmetric environment. In the present analysis, we describe the E1-E2 process 
by employing a cluster model of FeOe , where all the 3d and 4p orbitals of Fe atoms and the 2p 
orbitals of O atoms as well as the Coulomb and the spin-orbit interactions in the 3d orbitals 
are taken into account. We obtain an effective hybridization between the 4p and 3d states in 
addition to the ligand field on the 3d states through the hybridization with the O 2p states. 
On the basis of these frameworks, we clarify various symmetry relations to the E1-E2 process 
and the relation between the nonreciprocal directional dichroism and the anapole moment. 
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Furthermore, we numerically calculate the absorption spectra as a function of photon energy 
by diagonalizing the Hamiltonian matrix in the 3d^- and 4p3(i^-configurations. The shapes of 
magnetoelectric spectra as a function of photon energy are found similar to the experimental 
curves but their signs are opposite to the experiment in the pre-i^-edge region.^) The origin for 
the opposite sign is not known. Note that the magnetoelectric spectra in the optical absorption 
have been obtained in agreement with the experiment by using the same cluster model.®) 

This paper is organized as follows. In §2, we introduce a cluster model of FeOg. In §3, we 
describe the x-ray transition operators associated with Fe atoms. In §4, we derive the formulas 
of x-ray absorption, and present the calculated spectra in comparison with the experiment. 
The last section is devoted to concluding remarks. 

2. Electronic Structures 

2.1 FeOe cluster and 4p band 

In a FeOe cluster, we consider the Is, 3d and 4p states in the Fe atom, and the 2p states 
in O atoms. The Hamiltonian may be written as 



where 



+ + H^'-^'^ + H^'-^P, (2.1) 

ma 

^1/2 di/4 di/^ 

+ C3d XI ■ ^W<^')dlncjdm'a' 

mm' air' 

+ H,e- XCSU'dL^^ma', (2.2) 
maa' 

H^^ = Y.^'PjvaPJri'T, (2.3) 
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k»7'(T 

<b'^ = E<C'(^-Hl^'^w + H.c., (2.6) 
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^i.-4p ^ U^'-^^Y.P'iaP'nU'^'^'- (2-9) 

rfcrcr' 

The energy of Sd electrons can be described by H^*^ [eq. (2.2)] where dma represents an 
annihilation operator of a 3(i electron with spin a and orbital m (= —y^ —r"^ , yz, zx, xy). 
The symbol refers to the energy of 3d state with orbital m. The second and third terms 
of eq. (2.2) represent the intra-atomic Coulomb and spin-orbit interactions for 3d electrons, 
respectively. The matrix elements g {i'iV2'-, i^z^a) are expressed in terms of the Slater integrals 

F^, and [u stands for (m, cr)], and the spin-orbit coupling is Cad- We evaluate atomic 
values of F^, and Cad within the Hartree-Fock (HF) approximation,^) and multiply 0.8 to 
these atomic values in order to take account of the slight screening effect. On the other hand, 
we multiply 0.25 to the atomic value for since F° is known to be considerably screened by 
solid-state effects. The last term in eq. (2.2) describes the energy arising from the exchange 
interaction with neighboring Fe atoms, where (S)o-ct' represents the matrix element of the spin 
operator of 3d electrons. The exchange field Hxc has a dimension of energy, and is ~ k^Td^, 
with Tc ~ 250 K. Note that this term has merely a role of selecting the ground state by 
lifting the degeneracy and therefore the spectra depend little on its absolute value. The Hxc is 
directed to the negative direction of the c axis at Fel sites when the magnetic field is applied 
along the positive direction of the c axis. 

The energy of oxygen 2p electrons is given by H'^p where pjr^„ is the annihilation operator 
of the 2p state of energy with rj = x,y,z and spin a at the oxygen site j. The Coulomb 
interaction is neglected in oxygen 2p states. The H^y^^^ denotes the hybridization energy 
between the 3d and 2p states with the coupling constant t^'^^^. The energy of the 2p level 
relative to the 3d levels is determined from the charge-transfer energy A defined by A = 
E'i-EP + lbU{3d^) - 10?7(3d^) with E'^ being an average of E^. The multiplet-averaged d-d 
Coulomb interaction in the 3d^ and 3d^ configurations are referred to as U{3d^) and U{3d^), 
respectively, which are defined by [/ = F° - (2/63) - (2/63) F^. 

The H^'P represents the energy of the 4p states, where Pkrj'a annihilation operator 

of the 4p state with momentum k, rf = x, y, z, and spin a. The 4p states form an energy band 
e4p(k). The density of states (DOS) of the 4p band is inferred from the K-edge absorption 
spectra^) (see Fig. 2) . The H^y^^^ represents the hybridization between the 4p and oxygen 2p 
states with the coupling constant i^f The annihilation operator of the local 4p orbital p'^,^ 
may be expressed as p'^^,^ = (l/\/]Vo) Yl'kP'krj'a (-^o is the discretized number of k-points). 

The energy of the Is state is denoted as H^^ where represents the annihilation operator 
of the Is state with spin a. Finally, the interaction between the Is and 3d states and that 
between the Is and 4p states are denoted as H^^~^'^ and H^^~^'p, respectively. 

Table I lists the parameter values used in this paper, which are used in our previous paper 
analyzing the optical spectra in GaFeOa,^) and are consistent with the values in previous 
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Table I. Parameter values for a FcOg cluster in the 3d^ configuration, in units of eV. The Slater- 
Koster two-center integrals arc defined for the Fe atom at the center of the octahedron. 



F^{3d, 3d) 


6.39 




-1.9 


F^{3d, 3d) 


9.64 


{pd'K)2p^U 


0.82 


F^{3d, 3d) 


6.03 


{PP<y)2pAp 


3.5 


Cm 


0.059 


{PPT^)2p,4p 


-1.0 


A 


3.3 







calculations for Fe304.^°'^^^ 

2. 2 Ligand field and effective hybridization between the 4p and 3d states 

Fe atoms are assumed to be displaced from the center of the octahedron along the b axis. 
The shift S is 0.26A at Fel sites and — O.llA at Fe2 sites. We evaluate the hybridization 
matrices tmri'^^{j) and f^^^^iJ) for the Fe atom at the off-center positions by modifying the 
Slater-Koster two-center integrals for the Fe atom at the central position of the octahedron 
(Table I) with the assumption that {pda)2p,3d, {pdiT)2p,zd on d~^, and {ppcF)Ap,2p-, {ppT^)4p,2p on 
d~'^ for d being the Fe-0 distance. Thereby the ligand field Hamiltonian on the 3d states 
is given in second-order perturbation theory, 

^3.-3. ^ J2 &Ld^,^ + H.C., (2.10) 
mm' a 

with 

^it?'' = E C'^o-)4'-;^(j)/A, (2.11) 

where the sum over j is taken on neighboring O sites, and A = 3.3 eV is the charge transfer 
energy defined in §2.1. In addition to the ligand field corresponding to the cubic symmetry, 
we have a field proportional to (5^. The latter causes extra splittings of the 3d levels. Similarly, 
we can evaluate the effective hybridization between the 4p and 3d states in the form: 

H''-''=EK'rn''pU^'r + li.C. (2.12) 
rfmcr 

Here, the effective coupling is defined by 

iX'" = E - E'n, (2.13) 

where E'^p is the average of the 4p-band energy, which is estimated as E^p — E'^p 17 eV. 
The value of coupling coefficient i^f "^'^ is nearly proportional to the shift S of the Fe atom 
from the center of the octahedron. 
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2.3 Ground state 

Wc assume that Fe ions arc in the d^-configuration in the ground state, whieh will be 
denoted as \^g{d^)) with eigenenergy Eg(d^). We calculate the state by diagonalizing the 
Hamiltonian H-^^ + ^3d-3d^ where the exchange field Hxc and the displacement of Fe atoms 
from the center of the octahadron are different from Fel and Fe2 sites. The lowest energy 
state is characterized as ^Ai under the trigonal crystal field when the exchange field and the 
spin-orbit interaction are disregarded. The inclusion of these interactions could induce the 
orbital moment h{Lx), but its absolute value is given less than 0.004^. 

3. Absorption process on Fe 

The interaction between the electromagnetic wave and electrons is described by 



Hint = -^ J j(r)-A(r)d3r, 



(3.1) 



where c stands for the speed of light and j represents the current-density operator. The 
electromagnetic field A(r) for linear polarization is defined as 



A(r) 




27r/?2c2 



eCqe'^'' + H.c., 



(3.2) 



where Cq and Hujq are the annihilation operator and the energy of photon, respectively. The 
unit vector of polarization is described by e. We approximate this expression into a sum of 
the contributions from each Fe atom: 



H-mt = -- XI ^) ■ ^' ^) + 



(3.3) 



with 



j(q, i) 



A(q,e,z) 




''^■^'-'■')jnn'(r-r^)d3(r-r,) 



ecae 



The local current operator may be described by^^^ 



al(i)an'{i), 



jn„.(r-rj) = ^ [(V^;)^n' - 0;^V^!>„/] 

2 t 

- — A</.;0„, + — cV X [</.:s0„,]. 

mc mc 



(3.4) 
(3.5) 



(3.6) 



The integration in eq. (3.4) is carried out around site i, and a„(z) is the annihilation operator 
of electron with the local orbital expressed by the wave function ^^(r — r^). The charge and 
the mass of electron are denoted as e and m, and is the spin operator of electron. The 
second term in eq. (3.6), which describes the scattering of photon, will be neglected in the 
following. The approximation made by taking account of the process only on Fe atoms may 
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be justified at the corc-levcl spectra, since the core state is well localized at Fe sites. 

The absorption cxperiment^^ we analyse has been carried out on the geometry that the 
photon propagates along the a-axis with linear polarization, as illustrated in Fig. 1. Corre- 
sponding to this situation, it is convenient to rewrite the interaction between the matter and 
the photon in a form. 



i^int = E \/ E ^(q' ^' Ocqe^'i--^ + H.C., (3.7) 

q V i 

where the transition operator r(q, e, z) is defined as he • j(q, i)/e. The explicit expression of 
r(q, e, i) for the El and E2 transitions are given in the following subsections. 

3.1 El transition 

The transition operator r(q, e,i) for the El transition is obtained by putting e*^'^'""'"-'^ = 
1 in eq. (3.4). Therefore it is independent of the propagation direction of photon. For the 
polarization along the z-axis, the first term in eq. (3.6) is rewritten by employing the following 
relation 

./ '^n§;'^n'dh = -^(e„ - e„,) J (f)lz(l)r,'dh, (3.8) 

where e,„ and e„/ are energy eigenvalues corresponding to the cigcnstatcs 0„ and respec- 
tively. At the i^-edge, we assign the 4p states to (pn and Is state to (pn'- Hence the transition 
operator T^^ is expressed as 

r^^(q,e,z)=iS^^^iVf [p;t««.«-4»Pw«]- (3.9) 

icr 

where i runs over Fe sites. Non- vanishing values of the cocfRcicnts AT^^^'s arc given by Nrj = 
l/\/3 for the polarization along the i] (= x, y, z) axis, independent of the propagating direction 
of photon. The coefficient B^^ is defined by 



roo 

{e4p-eis) r^R4pir)Risir)dr, (3.10) 
Jo 

where Ruir) and R4p{r) are radial wave functions of the Is and 4p states with energy eis and 
e4p, respectively, in the Fe atom. The energy difference may be approximate as e^p — eis ~ 
hwq = hcq.^^^ Within the HF approximation B^^ is estimated as B^^ ^ 1.5 x 10~^cm • eV in 
the ls^3d^4p'^ ™^-configuration of an Fe atom.^) 

3.2 E2 transition 

The transition operator for the E2 transition is extracted from eq. (3.4) by retaining the 
second term in the ©XpclIlSlOIl G * 1 -|- iq • (r — Tj) -|~ • ■ • . Let the photon be propagating 

along the y-axis with the polarization parallel to the z-axis. Then we could derive a relation, 

/d TTh f' y ^ 
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+ 0;L,</>„,d3r, (3.11) 

where HLx is the orbital angular momentum operator. The last term should be moved into 
the terms of the Ml transition.^) At the K-edge, we assign the 3d states to ^„ and the Is 
state to (f)n' , respectively. Hence the transition operator T^^ may be expressed as 

T^\q,e,i) = -</i?^^2^iV^^2(q)[^l^^(^)^^(.) _^t(,)rf^^(i)]. (3.12) 
imcr 

When the photon is propagating along the y-axis, m is selectively yz with A^^^(q) = 1/(2a/T5) 
in the polarization along the z-axis, and m is selectively xy with N^y{c[} = l/(2\/l5) in the 
polarization along the x-axis, respectively. Note that a relation iV^^(— q) = —N^^{q) holds. 
The g^j.g defined by 

B^^ = (e3d - ei,) / r^R3dir)Ris{r)dr, (3.13) 
Jo 

where R^di^^) is radial wave function of the 3(i state with energy e^^d in the Fe atom. An 
evaluation within the HF approximation gives B^^ k, 1.7 x 10~^^cm^-eV.^) 

4. Absorption spectra 



Restricting the processes on Fe atoms, we sum up cross sections at Fe sites to obtain the 
absorption intensity /(a;q,e). Dividing it by the incident flux c/V, we have 

/(a;q,e) oc ^ ^ |(^^(i)|r(q, e, z)|*,(z))|2 

x5{hw^ + Ey- Ef), (4.1) 

where T(q,e,z) = T^^(q,e,i) +T^^(q,e,i), and g{{)) and |^'j(i)) represent the ground 
and the final states with energy Eg and Ef at site i, respectively. The sum over / is taken 
over all the excited state at Fe sites. 

In the Fe pre-K-cdgc region, the final states are constructed by perturbation theory start- 
ing from the states {\^rn(d^), ^sa )) in the configuration with the Is-core hole. Within the 
second-order perturbation of ij4p-3d^ they are given by 

|*/(0) = |$m(d^),ls£) 

+ ^|^n(d^),kr?V,lsa) _^ 

X ($n(d'),W(7,lsa|#^P-3'^|$^(/),ls£), (4.2) 
where Ef stands for the energy of the unperturbed state. It is defined as 

Ef = Em{d^) - eis - Ei^tils - d^), (4.3) 
where Ei^{,{ls — d^) is the interaction energy between the electron in the Is states and electrons 
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in the 3cf configuration. In the second term of eq. (4.2), E'^^^, is defined as 

Kw = ^"(^') - + ^4p(k) - -Eint(ls - 4pd^), (4.4) 

where £^int(ls — Apd^) is the interaction energy between the electron in the Is states and 
electrons in the 4p3d^ configuration. Symbols k.r)'a and Isa appeared in the bras and kets 
indicate the presence of an electron in the 4p state (kr/'cr) and the absence of a Is-core 
electron with spin a, respectively. Since the second term of eq. (4.2) is completely evaluated 
by \^mid^), Isa ) and Ef, the label / is specified by the m-th eigenstates of the configuration 
and the core-hole spin a. Notice that the lowest values of eqs. (4.3) and (4.4) correspond to 
the positions of the pre- and main-edges, respectively. 

The sum over k may be replaced by the integral with the help of the Ap DOS. In our 
numerical treatment, the position of the pre-edge energy is adjusted to the experimental 
value and the difference between the pre- and main-edges is chosen as the minimum of E'^^^, 
to be 12 eV higher than Eg{d^) — eis — Es^ai{ls — d^) . For simplicity, the explicit dependence on 
site i is omitted from the right hand side of eq. (4.2). From these wave-functions, we obtain 
the expression of transition amplitudes at site i by 

M(q, e, i; /) = M^^Q, e, i; /) + M^\c^, e, i; /), (4.5) 

with 

M^\q,e,z;/) = (*/(0|r^'(q, e, z)|^'g(i)) 

= ^($^(d6),lsa|#^P-3'^|cI>,(d^),kVa,lsa) 

rikr)' 

1 

^ 



id''), kr]'a, lsa\T^\ci, e, i) |$g (d^)) , (4.6) 



M^\ci,e,i;f) = {^f{i)\T^\ci,e,i)\^g{i)) 

= {<^mid''),lsa\T^\q,e,i)\<^g{d^)). (4.7) 
With these amplitudes, eq. (4.1) is rewritten as 

/(a;q,q,e) oc ^ ^ |M(q, e, i; /)p 



T/tt 



[hu;^ + Egid5)-Ef]^+T^' ^^-^^ 
where the (^-function is replaced by the Lorentzian function with the life-time broadening 
width of Is-core hole T = 0.8 eV. 

Now we examine the symmetry relation of the amplitudes. First, let the propagating 
direction of photon be reversed with keeping other conditions. In eq. (3.4), iq- (r — r^) is to be 
replaced by — iq-(r— rj). We know that N^^ has no dependence on q and that A'^,^^(— q) is equal 
to — A'^^(q). Since other conditions are the same, we have the new amplitudes (M^^)' = M^^ 
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and (M^2y ^ -M^"^. Second, let the local magnetic moment at each Fe atom be reversed 
with keeping the same shifts from the center of octahedron. The reverse of the local magnetic 
moment corresponds to taking the complex conjugate of the wave functions. Considering 
eq. (4.6) together with eq. (3.9), we have (M^^)' = — (Af^^)*. Similarly, considering eq. (4.7) 
together with eq. (3.12), we have {M^"^)' = (M^^)*. Third, let the shifts of Fe atoms from 
the center of octahedron be reversed with keeping the same local magnetic moment, which 
means the reversal of the direction of the local electric dipole moment. This operation causes 
the reversal of the sign of i^^p-sd^ However, no change is brought about to the 3d states in 
the 3(i^- and 3(i®-configurations, because the ligand field }{^<^^^<^ varies as 5^. As a result, we 
have the new amplitude (M-^^)' = -M^'^ from eq. (3.9) while (M^'^y ^ ^E2^ 

As already stated, the direction of the local magnetic moment could be reversed by revers- 
ing the direction of the applied magnetic field, since the actual material is a ferrimagnet with 
slightly deviating from a perfect antiferromagnet. Let /^(tjq, q, e) be the intensity for the ex- 
ternal magnetic field along the ibc axis. Then, from the second symmetry relation mentioned 
above, we have the average and the difference of the intensities as 

/(wq, q, e) = ^ [/+(u;q, q, e) + /-(wq, q, e)] 



oc 



J;; (E E [l^^'Hq, e, ^; /)P + |M^2(q, e, z; /)p 

r/TT 



[^q + Eg{d^) - Eff + r2 



and 

A/(wq, q, e) = /+(c^q, q, e) - /-(cjq, q, e) 



+ [M^2(q,e,i;/)]*M^i(q,e,i;/) 

r/TT 



(4.10) 



respectively. The M^^ and M^^ represent the amplitudes when the magnetic field is applied 
parallel to the c axis. For the average intensity [eq. (4.9)], the last term describes the if-edge 
spectra due to the E\ transition that the Is electron is excited to the 4p band. Although its 
main contribution is restricted in the main K-edge region, its tail spreads over the pre-K-edge 
spectra due to the life-time width. We see that the difference intensity [eq. (4.10)] is brought 



10/14 



J. Phys. Soc. Jpn. Full Paper 

about from the E1-E2 interference process. According to the symmetry relations mentioned 
above, it is expected to follow^) 

A/(a;q,q,e) cxq-^PiocW x MiocW, (4-11) 

i 

where Pioc(i) and Mioc(i) are the electric and the magnetic dipole moment of Fe atom at 
site i, respectively. Note that Pioc(i) is proportional to 5i[= (0,0,5)]. Then, the right hand 
side of eq. (4.11) is the sum of the local toroidal moment T{i) [= 61 x Mioc(?)]."'^^^ We have 
already derived the same form in the optical absorption spectra,*^) which is brought about by 
the El-Ml interference process. The spectra change their sign if one of the vectors among q, 
Ploc, or Mioc reverses its direction. 

Figure 2 shows the calculated average intensity I{ujq, q, e) as a function of photon energy 
huJci, in comparison with the experiment. The l,s-corc energy is adjusted such that the K- 
edge position corresponds to the experiment. The intensity at the i^-edge {huq > 7120 eV) 
mainly comes from the El-El process given by the last term of eq. (4.9). In the pre-i^- 
edge region (?2u;q ^ 7110 — 7115 eV), the tail of that intensity spreads due to the life-time 
broadening of the core level. In addition, wc have another contribution of El-El process 
through the term proportional to jM^-'^p, and that of the E2-E2 process through the term 
proportional to |M^^p. The latter is found larger than the former, giving rise to a small 
two-peak structure with a weak polarization dependence (the inset in Fig. 2). The El-El 
process through the term proportional to |M^^p is effective only on the noncentrosymmetric 
situation, because the 4p3(i^ls-configuration has to mix with the 3d^ls-configuration. Note 
that the El-El process could also gives rise to the intensity in the pre-iiT-edge region through 
the mixing of the Ap state with the 3d states at neighboring Fe sites. This process need not the 
noncentrosymmetric situation, and has nothing to do with the magnetoelectric spectra. It is 
known that a substantial intensity of the resonant x-ray scattering (RXS) spectra is brought 
about from this process at the pre-K-edge on LaMnOs,^^^ but the present analysis could not 
include this process because of the cluster size. 

Figure 3 shows the magnetoelectric spectra A/(a;q,q, e) as a function of photon energy 
hcoq. For polarization e parallel to the b axis, the calculated spectra form a positive sharp 
peak and then change into a negative double peak with increasing hujq. On the other hand, 
for polarization e parallel to the c axis, the calculated spectra form a negative sharp peak and 
then change into a positive sharp peak with increasing fkOq. The corresponding experimental 
curves look similar, but their signs are opposite to the calculated ones. We do not find the 
origin for the opposite sign. 

5. Concluding Remarks 

We have studied the magnetoelectric effects on the x-ray absorption spectra in a polar 
ferrimagnet GaFeOs. We have performed a microscopic calculation of the absorption spectra 
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Fig. 2. Average intensity 7(a;q,q, e) as a function of photon energy fiWq. Photons propagate along 
the positive direction of the a axis with polarization vector e along the b axis. The solid line repre- 
sents the calculated spectra. The broken line denotes the experimental data with the background 
intensity subtracted from the raw data given in rcf. [4]. The dotted line represents the 4pD0S; 
both the low-energy and high-energy sides are arbitrarily cut-off. 



using a cluster model of FeOe- The cluster consists of an octahedron of O atoms and an Fe atom 
displaced from the center of octahedron. We have disregarded additional small distortions of 
the octahedron. We have derived an effective hybridization between the 4p and 3d states as 
well as the ligand field on the 3d states by modifying the Fe-0 hybridizations due to the shifts 
of Fe atoms. This leads to the mixing of the 4p3(i^-configuration to the 3d^-configuration, 
and thereby to finite contributions of the E1-E2 interference process to the magnetoelectric 
spectra. We have derived the symmetry relations of the amplitudes M^^ and M^^, and have 
discussed the directional dichroism of the spectra. The cluster model used in the present paper 
is the same as the model used in the analysis of the optical absorption spectra in GeFe03 
and is similar to the model used in the analysis of RXS in Fe304.^^^ We have numerically 
calculated the magnetoelectric spectra as a function of photon energy in the pre-K-edge 
region. Although the spectral shapes are similar to the experimental curves, their signs are 
opposite to the experimental ones. The origin for the opposite sign has not been clarified yet. 
We would like to simply comment that the magnetoelectric spectra in the optical absorption 
are obtained in agreement with the experiment by using the same model.^^ 

The magnetoelectric effect on RXS has been studied experimentally^^^ and theoreti- 
cally^^' in GaFeOs. We think the approach used in the present paper is effective also 
to the analysis of RXS. Closely related to these studies, the magnetoelectric effects on RXS 
have also been measured, and theoretically analyzed^^^ in magnetite, where A sites are 
tetrahedrally surrounded by oxygens with the local inversion symmetry being broken. 
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Fig. 3. Difference of the absorption intensities A7(a;q, q, e) as a function of photon energy hujq when 
the magnetic field is apphed paraUel and antiparallel to the c axis. The sohd and broken hnes 
correspond to A/(a;q,q, e) where photons propagate along the positive direction of the a axis 
with polarization vector e along the b and c axes, respectively. Experimental data are taken from 
ref. [4] and denoted as filled (e || b) and open (e || c) circles, respectively. 
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